home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dsygs2.z / dsygs2
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDSSSSYYYYGGGGSSSS2222((((3333FFFF))))                                                          DDDDSSSSYYYYGGGGSSSS2222((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DSYGS2 - reduce a real symmetric-definite generalized eigenproblem to
  10.      standard form
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DSYGS2( ITYPE, UPLO, N, A, LDA, B, LDB, INFO )
  14.  
  15.          CHARACTER      UPLO
  16.  
  17.          INTEGER        INFO, ITYPE, LDA, LDB, N
  18.  
  19.          DOUBLE         PRECISION A( LDA, * ), B( LDB, * )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      DSYGS2 reduces a real symmetric-definite generalized eigenproblem to
  23.      standard form.
  24.  
  25.      If ITYPE = 1, the problem is A*x = lambda*B*x,
  26.      and A is overwritten by inv(U')*A*inv(U) or inv(L)*A*inv(L')
  27.  
  28.      If ITYPE = 2 or 3, the problem is A*B*x = lambda*x or
  29.      B*A*x = lambda*x, and A is overwritten by U*A*U` or L'*A*L.
  30.  
  31.      B must have been previously factorized as U'*U or L*L' by DPOTRF.
  32.  
  33.  
  34. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  35.      ITYPE   (input) INTEGER
  36.              = 1: compute inv(U')*A*inv(U) or inv(L)*A*inv(L');
  37.              = 2 or 3: compute U*A*U' or L'*A*L.
  38.  
  39.      UPLO    (input) CHARACTER
  40.              Specifies whether the upper or lower triangular part of the
  41.              symmetric matrix A is stored, and how B has been factorized.  =
  42.              'U':  Upper triangular
  43.              = 'L':  Lower triangular
  44.  
  45.      N       (input) INTEGER
  46.              The order of the matrices A and B.  N >= 0.
  47.  
  48.      A       (input/output) DOUBLE PRECISION array, dimension (LDA,N)
  49.              On entry, the symmetric matrix A.  If UPLO = 'U', the leading n
  50.              by n upper triangular part of A contains the upper triangular
  51.              part of the matrix A, and the strictly lower triangular part of A
  52.              is not referenced.  If UPLO = 'L', the leading n by n lower
  53.              triangular part of A contains the lower triangular part of the
  54.              matrix A, and the strictly upper triangular part of A is not
  55.              referenced.
  56.  
  57.              On exit, if INFO = 0, the transformed matrix, stored in the same
  58.              format as A.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDSSSSYYYYGGGGSSSS2222((((3333FFFF))))                                                          DDDDSSSSYYYYGGGGSSSS2222((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      LDA     (input) INTEGER
  75.              The leading dimension of the array A.  LDA >= max(1,N).
  76.  
  77.      B       (input) DOUBLE PRECISION array, dimension (LDB,N)
  78.              The triangular factor from the Cholesky factorization of B, as
  79.              returned by DPOTRF.
  80.  
  81.      LDB     (input) INTEGER
  82.              The leading dimension of the array B.  LDB >= max(1,N).
  83.  
  84.      INFO    (output) INTEGER
  85.              = 0:  successful exit.
  86.              < 0:  if INFO = -i, the i-th argument had an illegal value.
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.